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We give a new demonstration of the Magnetostatic Theorem of Parker (1994) within the 
model of Parker (1972) where the theorem was first proposed. This theorem shows how 
high electrical conductivity may result in significant resistive dissipation in a hydromagnetic 
plasma through the spontaneous formation of electric current sheets (Parker 1979, 1994). 
As a fundamental explanation of the ubiquitous presence of high-temperature, magnetized 
plasmas in the astrophysical universe, interest and debates over this theorem have appeared 



O \ in many publications over the years (Parker 1972, 1979, 1994, 2000, Rosner & Knobloch 1982, 



Tsinganos et al. 1984, van Ballegooijen 1985, 1988, Hahm & Kulsrud 1985, Moffatt 1985, 
Zweibel & Li 1987, Ng & Bhattacharjee 1998, Low 1990, 2006b, 2007, 2010, Bogoyavlenskij 
2000, Craig & Sneyd 2005, Low & Flyer 2007, Janse & Low 2009, 2010, Huang et al. 2009, 

CO 

Low & Janse 2009, Aly & amari 2010). The theorem has a clear physical statement but 

^p" 1 ! is difficult to express mathematically (Janse, Low & Parker 2010, hereafter referred to as 
O 

JLP). The formidable part of the theory is in the explicit description of the topology of 

cfl ! 

three dimensional (3D) solenoidal fields. We state the Magnetostatic Theorem in Section 



2, analyze in Section 3 some properties involved in boundary displacements of magnetic 
footpoints, and, finally, demonstrate the theorem in Section 4. We conclude the paper in 
Section 5. 

II. STATEMENT OF THE MAGNETOSTATIC THEOREM 



The induction equation 

describes the evolution of a magnetic field B in a perfectly conducting fluid moving with ve- 
locity v. Magnetic flux is frozen-into the fluid so that the magnetic flux surfaces (MFSs) and 
field lines (FLs) move as fluid surfaces and lines of fluid particles. Taking v to be continuous 
in space implies that these geometric surfaces and lines are continuously deformed. Topology 
as defined in standard mathematics textbooks deals with the properties of geometric objects 
that are invariant under completely general, continuous transformations, including transfor- 
mations that do not preserve any geometric metric. The fields at two instants of time are 
the topological transformation of one into the other with certain toplogical properties that 
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do not change. Ultimately, these properties are physically expressible as the conservation 
of magnetic flux across fluid surfaces (Low 2006a, LJ). Given any field B there is a set 
of topological properties we denote formally as T(B) that is invariant in time. Under the 
frozen-in condition, T(B) identifies a field physically for all time. 

Consider the infinite domain V : \z\ < L , L Q being a constant, in standard Cartesian 
coordinates. The boundaries z = ±Lo are perfectly-conducting rigid plates where the field 
has a fixed normal component, 

B z \ z=±Lo = F±{x,y). (2) 

The continuity of the tangential electric field at the plate, demanded by Maxwell's equations, 
imposes the boundary condition 

v|*=±i = 0, (3) 

assuming that F± ^ 0. We limit ourselves to fields that thread across V from one boundary 
plate to the other, that is, the case of subsystems of flux closed entirely within V is ex- 
cluded. Then the unchanging T(B) comprises the mapping of points on one plate to those 
on the other as pairs of boundary footpoints of FLs, the twist in all definable fluxtubes, and 
the twists among fluxtubes, topological features easily stated conceptually but difficult to 
describe in explicit mathematical terms (JLP). 

With a tenuous atmosphere like the solar corona in mind, we adopt the approximation 
that the non-magnetic forces of the fluid are negligible. Then the static equilibrium of the 
magnetized fluid is described by the force-free equations 

(V x B) x B = 0, (4) 
V B = 0, (5) 

stating that the field exerts no force on the fluid. Under the frozen-in condition, the field 
independent of its state has a distinct identity defined by T(B). So the construction of a 
force-free field in the above domain is not a classical boundary value problem, but a problem 
subject to both boundary condition (j2J) and a prescribed T(B). Any expression of this non- 
classical mathematical problem involves integral equations since topology is a global property 
(JLP). The Magnetostatic Theorem states that for most prescriptions of (F±(x,y),T(B)), 
equations (J4]) and (j3J) have no continuous solutions, which in classical analysis simply means 
no solutions. 
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The problem as posed is nevertheless physically meaningful, and it only means that its 
solutions must contain magnetic tangential discontinuities, the so-called weak solutions of 
mathematical analysis (Courant & Hilbert 1962, JLP). These are discontinuities of B across 
MFSs subject to continuity of B 2 everywhere, so that the integral version of equations (jl]) 
and §5§ is satisfied. The magnetic tangential discontinuities contain the current sheets and 
the continuity of B 2 ensures that these current-sheet surfaces are in macroscopic equilibrium. 

Let us concentrate on the case of a force-free field described by the infinite series 



where e is a constant small parameter, assuming that such a series is analytic with a finite 
radius of convergence independent of spatial position (Rosner & Knobloch 1982, JLP). 
We may think of some solution space of force-free states in which this series picks out a 
neighborhood around the uniform field of unit strength. Some of these states are produced 
from the uniform field under the frozen-in condition in the manner first considered by Parker 
(1972). Starting with the uniform field, first relax the rigid anchoring of the field at z — ±Lq, 
and then shuffle the magnetic footpoints by some continuous displacement to redistribute the 
boundary B z and introduce magnetic twist into the flux tubes in V . Now take the boundary 
plate as rigid again, freezing the new B z distribution and the new field topology created. 
The Magnetostatic Theorem then states that for most continuous footpoint displacements, 
the deformed field may not find a continuous force-free equilibrium. In the context of the 
infinite series ([6]), the claim is that there are footpoint displacements producing deformed 
fields whose force-free states are not found in the e-neighborhood of the uniform field. 

III. FIELD DEFORMATION BY BOUNDARY FOOTPOINT DISPLACEMENTS 

We take the footpoint displacements to be continuous. Any discontinuous displacement 
at the boundary naturally produces a current sheet extending from the boundary into the 
domain. Let us use the acronym CMFD for continuous magnetic footpoint displacements. 
We treat the properties of CMFDs, and relate them to a property of potential fields studied 
by Low (2007) and Janse & Low (2009), hereafter referred to as L07 and JL09. This prepares 
us for the demonstration in Section 4. 



B(x,y,z) = z + b(x,y,z) 



oo 




(6) 



n=l 
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A. Field Topology and Footpoint Motions 



Take the z component of the induction equation (CQ) for a velocity v = u = (u x .u y , 0) on 
some plane of constant z, 

Ft R 

^ + Vj. • (S.u) = 0, (7) 

where the subscript _L indicates partial differentiations in the constant z plane. If u is 
given, equation (J7J) poses an initial value problem for B z (x, y, z, t) evolving from some initial 
distribution B z (x,y,z, 0) at t = 0. If we know B z (x,y, z,t), instead, equation ([7]) poses a 
static problem at each instant of time to determine the velocity u. 
Express the velocity in the general form 

u = V ± (p + V ± x (ifjz) , (8) 

separating the compressible irrotational from the incompressible rotational parts, described 
by and if), respectively. If we set if) = 0, the fluid merely compresses the magnetic flux 
threading across the constant-z plane without twisting the field. Similarly, if we set <f> = 0, 
the fluid twist the field threading across that plane without compression. In the latter, B z 
can still change with time because of the frozen-in flux transport in the constant- z plane. 
To illustrate, set • u = to obtain 

dB 

+ u ■ V ± B Z = 0, (9) 

offering the possibility of twisting the field without changing the ^-distribution by an in- 
compressible rotational displacement with if)(x, y, t) being a strict function of the unchanging 
B z . 

Consider a given B z (x, y, z, t), with B z (x, y, z, 0) = H (x, y) at t = and B z (x, y, z, ti) = 
Hi(x,y) at t = ti, at two chosen times, keeping in mind we are dealing with events in a 
constant-z plane. That evolution can be accounted for by a purely compressive irrotational 
flow. Set if) = and then solve 

rlB 

^ i + V ± -( J B 2 V ± 0) = O, (10) 

as an elliptic partial differential equation for (f>. The solution is unique under suitable 
physically reasonable boundary conditions for the unbounded constant-^ plane. In other 
words, the evolution from B z = H at t = to B z = Hi at t = t\ can be achieved through a 



purely irrotational footpoint displacement. No magnetic twist is introduced in this process. 
Upon arrival at the state B z — H\ at t — t\, magnetic twist can be put into that field 
subsequently, if so desired, without changing B z beyond this point in time. This is done by 
an incompressible rotational displacement x/j = ip(Hi,i) that leaves B z unchanged in time. 

The details of B z (x,y,z,t) in the period < t < t\ determines the precise form of the 
associated velocity u = V±(p in the above process. There is an infinity of prescribable 
B z (x, y, z, t) sharing a pair of initial and final states (Hq, H\). Hence, there is an infinity of 
compressible, irrotational CMFDs that can take B z from initial to final prescribed distrib- 
tions. Now let us distinguish between the Eulerian velocity of the fluid from the Lagrangian 
velocity of a fluid particle. The CMFDs are the Lagrangian displacements of actual fluid 
particles on the constant-z plane. Different Eulerian velocities may produce the same fi- 
nal ^-distribution, but their respective actual CMFDs are quite distinct. In other words, 
the footpoints on the constant-;? plane are displaced continuously in different manners al- 
though all the different Eulerian velocities bring about the same final .Redistribution. It 
should be emphasized that this is possible only if the footpoints are free to move in the full 
two- dimensionality of the constant- z plane. 

Returning to the Parker two-plate domain with an initial uniform field, we draw two 
important conclusions. Firstly, if the CMFDs at the boundaries z = ±Lo are incompressible, 
then it is obvious that the distribution B z = 1 remains at those boundaries. The field merely 
acquires a twist. The different terms in the series (|6]) must then satisfy b njZ = at z = ±Lo 
for all n; see the study Low (2010). 

If the CMFDs at the boundaries z = ±Lo are compressible but irrotational, then B z must 
change at those boundaries. Thus, we pose the boundary conditions 

b n ,z\z=±L = fn,±(x,y), (H) 

for n = 1,2, 3,..., where f Ut ± are prescribed. In other words, we have the expansion 

oo 

^±(x,?/) = l + ^e n /n, ± (x,l/), (12) 

n=l 

to describe the final B z at the two boundaries produced by an imposed CMFD. We had 
simply treated e as a constant parameter until now. The free specification of B z at the 
boundaries may be taken to be the physical origin of that parameter. 
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When we displace the footpoints of an initial uniform field at both plates with a rotational 
motion, i.e., ip ^ 0, it does not mean that we would have twisted the field. It is the relative 
motion of the pair of footpoints of a FL that determines whether a twist has been built into 
the field. To avoid this qualification, we shall henceforth consider deforming the uniform 
field with the imposed CMFDs taken only on z = — L , that is, the magnetic footpoints are 
not displaced on z = L where there is no change in the ^-distribution. There is no loss of 
generality to the principal physical point we wish to make. In this case, a rotational CMFD 
twists the uniform field whereas a compressive, irrotational (ip = 0, <j> ^ 0) CMFD introduces 
no twist on all scales. What is meant by untwisted is in this case unambiguously defined 
for the field deformed with the latter CMFD. Moreover, there is an infinity of compressive, 
irrotational CMFDs that can take the uniform field to deformed states, all sharing the same 
^-distributions at the two boundaries but not topologically equivalent by virtue of their 
distinct footpoint connectivities. 

B. The Janse-Low Result 

The possibility of an infinity of untwisted fields in some given domain, all fields sharing 
the same boundary flux distribution but having different topologies, was encountered in 
studies L07 and JL, treating the domain between two concentric spheres and the upright 
cylindrical domain, respectively. Let us briefly summarize the JL result. 

Consider a potential field B pot (L) in form of a unidirectional flux entering an upright 
cylinder of length L and radius Rq, from one cylinder end z = — L and exiting at the 
other end z = L. This potential field is tangential at R — Rq, using the usual cylindrical 
coordinates. The boundary flux distributions B z , positive definite in this case, at z — ±L 
define this potential field uniquely. By integrating 



with B = B pot we obtain each of its FLs described by x(z) and y(z), two coordinates of a 
point along the FL in terms of the third z as the independent variable. These FLs define 



dx 
dz 

dy 

dz 



B x (x,y,z) 
B z (x,y,zY 
B y (x,y,z) 
B z {x,y,zY 



(13) 
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the map M [B pot (L)] : (x B , Ub, -L) ->■ (x T , Vt, -L): 

x T = x B + [ Z=L B A<z),y{z),A dz 

B J z =-L B z [x(z),y(z),z] 
f z=L B v [x(z),y(z),z] , 

with B = B pot , relating a footpoint (xb,Ub, —L) to the footpoint (xt,Ut, —L). It was 
discovered in JL that the map M. [B pot (L)] is generally dependent on L, the exceptions 
being cases containing special symmetries. These symmetries include the obvious case of 
axisymmetry as well as field that are symmetric about the z = plane, for examples. This 
Independence of the footpoint map is a 3D effect. 

Consider deforming continuously a potential field B poi (Li), along with its domain 
L = L\ 7^ Lq, under the frozen-in condition into a non-potential field B^(Li;Lo) to fit 
into the domain L = L , treating the cylinder ends as rigid conductor. Throughout our 
discussion we assume no change in the cylinder radius. In this deformation the footpoint 
map M is a topological invariant. There is no way for the deformed field 'B d f(L 1 ; L ) to 
assume the potential state B po t(L ) by some suitable further deformation without also mov- 
ing the footpoints at the boundary, because their footpoint maps M. are not the same. 
With its invariant topology, the deformed field must find some force-free state distinct from 
B pot (L ) and the arguments in JL concluded that the force-free state must contain tangential 
discontinuities. 

For our deductive analysis here, it is possible to be rigorously specific about a concept 
of an untwisted field without having to describe its topological properties explicitly. We 
call all potential fields, unique in their respective (simply-connected) L-domains, untwisted; 
see JLP for an important discussion of multiply-connected domains. This is based on the 
fact that the magnetic circulation vanishes for the potential field along any closed curve 
in the simply-connected domain, a well known property. For the field B pot (L) = z + b 
with |b| << 1, every flux tube extending from z = —L to z = L shows no circulation 
around the approximately straight tube with an axial flux conserved along it. If this field 
is deformed and kept within the regime |b| << 1, magnetic circulations, quantities that are 
not invariant under the frozen-in condition, can be created in opposite signs along the flux 
tube. Conversely, the opposite circulations found at any one time along a flux tube can be 
removed by a suitable mutual cancellation of all of them along the tube. One thing that 
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is forbidden is to deform any of these potential fields to end up with a positive or negative 
definite along the tube. 

With this definition of untwistedness, a specific L = L domain contains an infinity 
of untwisted fields sharing the same boundary flux distribution, namely, fields that are 
deformed from the potential fields B pot (L) originally occupying domains with L ^ L but 
deformed to fit into the L = Lq domain. The claim is these deformed fields must find 
force-free states containing tangential discontinuities. 

The deformation of the uniform field in a L ^ L domain into a potential field of the form 
B pot (L) = z + b does not involve rotational footpoint displacements because the magnetic 
circulation of the potential field vanishes for all closed curves. The same particular irrota- 
tional CMFD can also be applied to the boundaries z = — L of the L domain containing a 
uniform field. The deformed field so produced is topologically identical to B po t(L), that is, 
they have the same footpoint map. This clarification assures us that the different maps M. 
of all the different untwisted fields in the Lq domain are produced by irrotational CMFDs. 

IV. THE e-NEIGHBORHOOD OF THE UNIFORM FIELD 

We first solve the magnetostatic equations for the field given by e-series (JSJ), without 
consideration of field topology, and then examine the implications of such analytic solutions 
in terms of the field topologies they exhibit. 

A. Perturbational Analysis 

Rewrite the force-free equation (J5J) in the form 



2(B- V)B- VB 2 = 0, 



(15) 



and substitute for the field given by series (jSJ) to first and second orders, 

d 

- V6i, 2 = 0, 

oz 

2 (bi • V) bx + 2^-b 2 - V (6? + 2b 2 , z ) = 0, 



(16) 



(17) 



subject to the solenoidal conditions 



V • bi = V • b 2 = 0, 



(18) 
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and boundary conditions (fTTj) on b\ )Z and 62,2- No other boundary conditions are needed in 
this consideration that has not yet introduced the question of field topology. 
The first order equation (TIB]) subject to the solenoidal condition implies 



vV 



0, 



(19) 



which has a unique solution satisfying the prescribed boundary values at z — ±L . Here 
and elsewhere we assume the far boundary condition that all departures from the uniform 
field vanish at x 2 + y 2 — > 00. The solenoidal condition shows that the first order magnetic 
field has the form 

bi = V$ 1 (x,y i z) + a 1 (x,y), (20) 

where ai(x, y) = (ai )X , a>i,y, 0) is an arbitrary solenoidal vector and V$i is the unique poten- 
tial field whose z component accounts for b\ jZ computed above. If ai = 0, bi is an untwisted 
potential field. 

First note that the z component of equation ( 1TT1) is an identity whereas the other two 
components take the forms 



db 



2,x 



db 



'2.z 



dz 
db 2 , y 
dz 



dx 
db 2>z 
dy 



— u 



— u 



9$! 

dx 
d<h_ 
dy 



■v 



0, 
0, 



(21) 



where we have expressed hi using equation ( 120]) and introduced u = ^-^r-~^f-- Eliminating 
62,2 and &2,2/ by using the solenoidal condition in favor of 62,2 gives the Poisson equation 



V 2 6 2 ,* + — 



d_ 

dx 



w I -7: h a 1>x 

dx 



d_ 

dy 



d<f>t 
- I -Q- + ai, y 



0. 



(22) 



The solution b 2;Z subject to its boundary values at z — ±L is unique. This solution, 
substituted into equations (12T1) . determines the ^-dependent parts of b 2 , x and b 2y y uniquely. 
We are free to linearly superpose this second order magnetic field, denoted by b^x, z), 
with an arbitrary solenoidal vector a^ix, y) = (^2,^,02.^,0), and still have a valid solution. 
So we express the solution as 



b 2 = bl(x,y,z) +a 2 (x,y). 



(23) 



The free vectors a x and a 2 describe the twist in the net field ([6]) up to second order. 
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For example, if ai ^ 0, this solenoidal field by itself has FLs that are the same closed 
curves on all constant-z planes. This property follows from the boundary condition that 
the perturbations all vanish at x 2 + y 2 — > oo. The implied magnetic circulations around the 
principal flux in the z direction of any flux tube generally cannot mutually cancel. Thus, to 
first order in e the FLs wind with a fixed sense about flux tubes directed in the z direction. 
An untwisted field corresponds, at first order in e, to ai = 0, in which case bi = V$i, 
a pure potential field. For such a field, co = 0, and 62^ is potential by equation ([22]) . so 
that hl(x,y,z) is a potential field. In this case, the free, z- independent vector a 2 imposes 
a magnetic twist at the second order in e. Thus, if there is no twist, we must set a 2 = 0. 
The magnetic field B = 1 + ehi + e 2 b 2 is thus potential to second order. By induction, the 
perturbations to all orders are all potential and so the net field is potential. 

B. Footpoint Displacements 

Substituting a e-series solution into equation (TTJJ defines the footpoint map Ai of the 
points on z = —L Q to the points on z = L . Expanding in e, we obtain the formal expressions: 

00 

x T = x B + ^2e n Ax n (x Bl y B ), 

71=1 
OO 

Vt = ys + J^VA^XB,^), (24) 

n=l 

where Ax n and Ay n are complicated functions defined through the definite integrals with 
respect to z, straightforward to work out from equation f[T4"|) . The topological magnetostatic 
problem of Parker (1972, 1994) is posed by giving the footpoint map Ai, among other speci- 
fications of the topology of the force-free field to be constructed. In the e-series formulation, 
this means prescribing the functional forms of Ax n (xB,yB) and Ay n (xB,yB), and seeking 
the fields h n . This is obviously a generally intractable problem. Fortunately, we can derive 
a clear conclusion by the following particular problems. 

Suppose, we take a potential field in V of the form B = z + eb pot , where b pot is a given 
potential field independent of e taken to be small so that |eb pot | << 1. What we have is a 
series solution with bi = b pot , and b n = 0, for n > 2. In other words, the series ([6]) has 
only a finite number of terms, just 2 terms. Equations (|24l) giving the footpoint map Ai 
remain a pair of infinite series. Given h pot , each terms of these two infinite series can be 
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computed explicitly. Turning the problem around, specifying this footpoint map M. to pose 
the topological problem leads to its solution B = z + eh pot . This is a case of a deformed 
field whose prescribed field topology allows it to attain a force-free state, namely, h pot . 

We can also carry out this computation for a footpoint map Ai belonging to a deformed 
field B = z + eb^, not in equilibrium, obtained as followings. Construct first a potential 
field hpot^L in a domain V(L) : \z\ < L with L ^ L , subject to -Redistributions on z = ±L 
that are identical to the Redistributions of b pot on z = ±L . Compute the footpoint- 
map functions Ax n (xB,UB) an d Ay n (xB,UB)- Now continuously deform z + b pot)L along 
with the domain V(L) : \z\ < L, to obtain b d f that fits into V(L ). The footpoint-map 
functions Ax^xbiUb) an d Ay n (xB,VB) remain unchanged and the field remains untwisted. 
In demanding that this deformed field be further deformed under the frozen-in condition to 
be force free, we conclude from our perturbational analysis that it cannot find a continuous 
solution in the e-neighborhood of the uniform field. If it could, that force-free state has to 
be the unique potential field h pot . That is not possible because the two fields have different 
footpoint maps. This demonstrates the Parker Magnetostatic Theorem. 



V. CONCLUSION 



We have returned to the model of Parker (1972) to demonstrate the Magnetostatic The- 
orem described in generality in Parker (1994). The general theory is concerned with the 
infinite set of field topologies 7/// of all continuous force- free fields admissible in a physical 
system. This set is generally a true subset of the set T a u of the topologies of non-force-free 
continuous fields admissible in the same system, that is, 7/// C T a u- Thus a non-force-free, 
continuous field with a topology not found in 7/// must develop tangential discontinuities as 
the only way to become force-free as it perserves its topology under the frozen-in condition. 

We have demonstrated the Magnetostatic Theorem by treating the untwisted fields ad- 
missible in a physical system. The problem is simpler than the one posed by twisted fields. 
In 3D systems, the set T^ twd of topologies of untwisted non-force-free fields is generally 
infinite. In contrast, a 2D field with no neutral points; for example, an axisymmetric field, 
has a unique footpoint map, see JL and JLP. There is only one force-free state for an 
untwisted field in the 3D system, the unique potential field. Hence the set of untwisted 
topologies Tf L f l c wd realizible in a continuous force-free state is a set with a single member, 
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and, obviously, 77/} C T^f , thus demonstrating the theorem. 

A similar demonstration was made in the studies L07 and JL, but those studies depend on 
arguing around topological ideas that are extremely difficult to pin down mathematically. In 
the present paper, the study of the force-free fields in the e-neighborhood of the uniform field, 
treated by Parker (1972), simplifies the general problem. It is the ordering of the terms by 
magnitude in the series expansion, beginning with the dominant uniform field, that allows a 
clear sorting out between what is an untwisted field from those that are twisted. Basic to this 
analysis is the property that different footpoint connectivities can be created by footpoint 
displacements imposed on a uniform field, without twisting the field but all leading to the 
same boundary flux distribution. 

Counter arguments about the result in JL, like those recently presented by Aly Sz Amari 
(2010), need to be reexamined in light of this new result. We mention the following open 
questions that qualify the new result. The fact that we could not find continuous solution 
in the e neighborhood of the uniform field does not preclude the existence of the desired 
continuous solutions represented by some other small-parameter expansions; see the inter- 
esting development of van Ballegooijen (1985) and the discussions in Low (1990, 2010), 
Parker (1994), and JLP. Another possibility is that, despite the small-parameter nature of 
the CMFDs expressed by equation (1241) . continuous force-free solutions may exist in V that 
do match the prescribed footpoint displacements but can only be attained via deformations 
of finite amplitudes inside the domain. These questions should not detract from the cen- 
tral feature of our result, that the continuous force-free fields in the e-neighborhood of the 
uniform field are topologically of a more restricted variety than the non-force-free fields. 
The Parker Magnetostatic Theorem makes the fundamental point that this hydromagnetic 
feature is general. 

The National Center for Atmospheric Research is sponsored by the National Science 
Foundation. 
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